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Abstract 

Anomalous coarsening in far-from equilibrium one-dimensional systems is 
investigated by applying simulation and analytic techniques to minimal hard 
core particle (exclusion) models. They contain mechanisms of aggregated 
particle diffusion, with rates e<l, particle deposition into cluster gaps, but 
suppressed for the smallest gaps, and breakup of clusters which are adjacent 
to large gaps. Cluster breakup rates vary with the cluster length x as kx a . 
The domain growth law (x) ~ (et) z , with z = 1/(2 + a) for a > 0, is ex- 
plained by a simple scaling picture involving the time for two particles to 
coalesce and a new particle be deposited. The density of double vacancies, at 
which deposition and cluster breakup are allowed, scale as 1/ [i(et) 2 ]. Numer- 
ical simulations for several values of a and e confirm these results. A fuller 
approach is presented which employs a mapping of cluster configurations to a 
column picture and an approximate factorization of the cluster configuration 
probability within the resulting master equation. The equation for a one- 
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variable scaling function explains the above average cluster length scaling. 
The probability distributions of cluster lengths x scale as P(x) = j^g{y), 
with y = x/(et) z , which is confirmed by simulation. However, those distri- 
butions show a universal tail with the form g(y) ~ exp ^— y 3//2 ^ , which is 
explained by the connection of the vacancy dynamics with the problem of 
particle trapping in an infinite sea of traps. The high correlations of sur- 
viving particle displacement in the latter problem explains the failure of the 

independent cluster approximation to represent those rare events. 
PACS numbers: 05.50.+q, 05.40.-a, 68.43.Jk 
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I. INTRODUCTION 



Domain growth in far from equilibrium systems is a subject of increasing interest due to 
the large number of applications, such as phase separation of mixtures, dynamics of glasses 
and island coarsening after thin film deposition [2,3,1]. In these systems, their dynamics is 
responsible for bringing them to steady states, while external agents act to drive them out of 
equilibrium. Many statistical models exhibit simple, universal domain growth laws, which 
are found in some real systems, but there is much interest in models with slow coarsening and 
with continuously varying growth exponents, for instance due to their potential applications 
to glassy systems [2-4]. 

In this paper, we will consider one-dimensional models with particle deposition and 
diffusion, reversible aggregation to clusters and mechanisms of cluster breaking which show 
such a variety of domain growth laws. Cluster breaking may be an effect of internal stress 
and was previously considered in studies of island growth in sub monolayers [5]. In a real 
system subject to external pressure but with some type of geometrical frustration other 
than those observed in island growth, it is expected that cluster breaking will compete with 
mechanisms of densification, these ones to be represented by vacancy filling (deposition of 
new particles). However, the onset of those processes depends on the formation of large 
vacancies due to the (slow) diffusion of aggregated particles. The coarsening exponents of 
those systems will be shown to vary with the exponents in the scaling of the cluster break 
probabilities, although the cluster lengths distributions are universal. Consequently, these 
one-dimensional statistical models, although not related to a specific real problem, reveal 
some interesting features that may help to understand complex three-dimensional systems, 
with the advantage of being more tractable both analytically and numerically. 

The models presented here are nontrivial extensions of those analyzed in a recent pa- 
per [6], which include particle diffusion, reversible aggregation to clusters and deposition 
mechanisms. In the original model, hard-core particles in a one-dimensional lattice have 
diffusion rates r = d when they were free (i. e., they have two empty nearest-neighbor sites) 
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and r = e ~ e~ E l T when they have one occupied nearest neighbor site, with r <C d (Fig. la). 
The deposition rate is F = 1, in units of monolayers per time step, and is restricted to sites 
with at least one empty nearest neighbor (Fig. lb), i. e. a site of a double or larger vacancy. 
These dynamical rules were motivated by the Clarke- Vvedensky model and related models 
of thin films or submonolayer growth [7] , but included effects of geometrical frustration that 
forbid filling of single vacancies. Domain growth in the form (d) ~ e - 1 / 2 t 1 / 2 was predicted 
analytically and from simulation [6]. The same model without deposition and in the limit 
f < 1 showed domain growth as t 1 ^ before approaching a steady state [6]. 

Here, in addition to the processes of the original model (Figs, la and lb), we will consider 
the competition between deposition of new particles and the breaking of a neighboring cluster 
when a double (or larger) vacancy appears. In this model, a cluster of length x with two 
vacancies at one of its sides (where deposition may also occur) may break in two pieces, at a 
random internal position, with rate kx a , where the exponent a is a tunable parameter and 
A; is a constant amplitude. The process is illustrated in Fig. lc for a cluster with x — 6, with 
a total of 5 internal points for its separation into two pieces. Our focus is the non-trivial 
case a > 0, for which a-dependent coarsening exponents are obtained. 

Notice that this model is significantly different from other models which involve coagula- 
tion or breakup of clusters with size-dependent rates with variable exponents [8-11]. Here, 
the clusters slowly gain mass by (non-biased) diffusion at the expense of neighboring clusters, 
while breakup is subject to the availability of free space for its expansion. 

Scaling laws for the average cluster length (x) in the form 

(x) ~ f (1) 

will be obtained, where z is called the coarsening exponent. In the models with the cluster 
breaking mechanism, the exponent z can be continuously tuned from to 1/2 by varying 
the scaling exponent a of the rate of cluster breaking. This result is predicted by a scaling 
theory that also described the cluster growth laws of Ref. [6] and is confirmed by numerical 
simulations with very good accuracy. It is also possible to describe such systems in terms 

4 



of interval probabilities, which allows the analytic investigation based on an independent 
cluster approximation. Using this method, we also predict the coarsening exponents of the 
model. 

There is also much interest in knowing the distributions of cluster lengths in such anoma- 
lously growing systems. Those distributions are also calculated numerically and show an 
universal (a-independent) form P (x) ~ exp x 3 ^ 2 ^j , despite the fact that the coarsening 
exponents do depend on a. The analytical treatment of the model in the independent cluster 
approximation gives instead an a-dependent power in the exponent (— x 1+a ^ 2 ). However, 
it can be shown that the dynamics of large clusters is related to the problem of particle 
diffusion in an infinite sea of mobile traps [12-14], where the holes between clusters in the 
former problem correspond to the particles in the latter. This connection leads to the above 
universal cluster length distribution and explains the failure of the independent cluster ap- 
proximation to predict those distributions. 

This paper is organized as follows. In Sec. II we present the scaling theory and obtain 
the coarsening exponents. In Sec. Ill we present the results of simulations for the time 
dependence of average cluster lengths and density of double vacancies. In Sec. IV, we 
map the cluster configurations to a column picture and determine the master equation 
with the independent interval approximation to the joint cluster length probability. In Sec. 
V, we discuss the cluster length probabilities, comparing numerical results, the analytical 
prediction of the independent interval approximation and the connection to the problem of 
one diffusing particle in a sea of moving traps. In Sec. VI, we summarize our results and 
present our conclusions. 

II. MODELS AND SCALING THEORIES 

Here we will define our models and estimate coarsening exponents using scaling argu- 
ments along the lines of Ref. [6], which were previously adopted in the analysis of related 
systems by Evans [3] and introduced in the analysis of domain growth in magnetic systems 
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by Lai et al [15] and Shore et al [16]. 

First we consider the original model of deposition and diffusion presented in Ref. [6] 
(Figs, la and lb). 

For simplicity, we will refer to the average cluster length as x. In Fig. 2a, we show a 
configuration with clusters of lengths typically of order x, named A, B and C, with single 
empty sites (single vacancies) between them. Deposition is not allowed at those vacancies, 
as well as at the other single holes separating A and C from other neighboring clusters. 
Suppose that the length of A tends to increase in time, while the length of B decreases, as 
shown in Fig. 2b. This evolution is equivalent to the diffusion of the vacancy between A 
and B, which gets closer to the vacancy between B and C (diffusion of this vacancy was not 
shown only to simplify the illustration). Finally, the length of A will increase by x when A 
and B coalesce, as shown in Fig. 2c. At that time, a new particle is deposited in an empty 
site of the double vacancy, which is also shown in Fig. 2c. 

The deposition process occurred after all particles of cluster B have detached from it 
and aggregated to cluster A. This is equivalent to displacement of order x of the diffusing 
vacancy between A and B. Notice that in a configuration with two neighboring empty sites 
(Fig. 2c), the probability of deposition (fixed rate F — 1) is much larger than the probability 
of diffusion of an aggregated particle (e <C 1), thus it is highly improbable that a diffusion 
process will follow the formation of a double vacancy. Consequently, the time necessary for 
coarsening of two clusters of length x is of the order of the time for a diffusing vacancy to 
move a distance x, which is At ~ e~ 1 x 2 . In the configuration after deposition, the average 
cluster length is increased by x, since it was equivalent to the merging of cluster B into 
cluster A. Consequently, the average cluster length increases as 

dx Ax x e 
~dt ~ At ~ e~ l x 2 ~ x' 

Integrating Eq. 2, we obtain x ~ e^t 1 / 2 , in agreement with the analytical results of the 
independent cluster approximation and simulation data [6]. 

Now consider the problem with the cluster breaking mechanism. It is assumed that a 
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(2) 



cluster of length x breaks with rate 

r cb = kx a (3) 

only when there is enough space available, i. e., when there is more than one empty site 
at one of the sides of the cluster. This process is illustrated in Fig. lc, in which a cluster 
of length x = 6 may break at five different internal points, which gives five possible final 
configurations if the breaking occurs - three of them were shown in Fig. lc. 

Following the argument of the original model, after the diffusion processes leading to 
the onset of a double vacancy (Fig. 2c), there are three possibilities: a) a new particle is 
deposited with rate F — 1; b) the cluster at the right side of the double vacancy breaks; c) 
the cluster at the left side breaks. Since the two neighboring clusters have lengths of order 
x, deposition will occur with probability Pd ep ~ l/x a for a > and i C 1. Otherwise, one 
of the clusters adjacent to the double vacancy will break in two pieces. 

If a new particle is deposited, then the average cluster length will increase by x, such as 
in model I. On the other hand, in the most probable case of an adjacent cluster to break, the 
average cluster length will not increase. Consequently, the average cluster length increases 
as 



dx Ax 1 x e 

~di ~ dep At ~ ^F 1 ^ = ( ' 



Integrating Eq. 4, we obtain 



x ~ et z , z = . (5) 

2 + a y ' 



For a < 0, the above reasoning leads to Pd ep ~ 1 for large x and, consequently, z— 1/2, 
as in the original model. 

For a > 0, Eq. 5 shows that the coarsening exponent can be continuously tuned from 
z = to z = 0.5 by changing the cluster breaking exponent. 

The same arguments can be used to predict the density of double vacancies, poo (t), a 
quantity which plays a important role in the analytic calculations of Sec. IV. The total 



rate at which deposition and cluster-breaking occur after the formation of a double vacancy 
(Fig. 2c) is of order x a , for a > 0. Consequently, that vacancy will survive during a time 
l/x a , while it takes a time of order At ~ e^x 2 to be formed (Figs. 2a-2c). Consequently, 
one double vacancy between two clusters of length x typically survives during a fraction 
(l/x a )/(e~ 1 x 2 ) of the total time. To obtain its density we must divide this quantity by the 
average cluster length x, which gives 

(l/x a ) 1 e et~ z 
Notice that this time decay does not involve the diffusive factor et alone. 



III. SIMULATION RESULTS FOR AVERAGE LENGTHS AND DENSITIES 

We simulated the model for cluster-breaking exponents a = 0.5, a — 1, a — 2 and a — 4, 
with e = 0.001 and, in some cases, also with e = 0.01. Deposition rates, free particle diffusion 
rates and amplitudes of cluster breaking rates were F — 1, d — 1 and k — 1, respectively, in 
all simulations. Lattices lengths were L = 5 x 10 4 . The maximum simulation times were 10 6 
for the smaller value of a and 10 8 for the largest one. Thus, at all times, the whole lattice 
contained more than 500 clusters, which was essential to avoid finite-size effects. Indeed, 
several results were compared with those in lattices of lengths L = 10 4 and confirmed the 
absence of significant finite-size effects. The number of different realizations varied from 
5000 for the smaller a to 500 for the largest one. 

The scaling of the average cluster length (x) on e and t is confirmed in Fig. 3 for a = 1, 
by showing (x)/ (et) 1 ^ 3 for e = 0.01 and e = 0.001 as a function of 1/ (et) 1 / 3 . The data collapse 
confirms the expected scaling on e, the convergence to a finite value as t — > oo confirms the 
scaling on t and the variable in the abscissa suggests a constant correction term to Eq. (4). 

In Fig. 4a and 4b we show (x)/(et) z versus l/(et) z for a = 0.5 and a — 2, respectively, 
with the values of z given by Eq. (4) and e = 0.001. These results confirm the validity of 
the scaling theory of Sec. II for the average cluster length. 
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Since double vacancies are very rare, the accuracy of their density is much lower than the 
accuracy of (x) . Consequently, averages over relatively large time regions were necessary to 
study the long time evolution of that density. The average values agree with the predicted 
scaling (Eq. 6), as shown in Fig. 5 for a = 1 (with two values of e). Notice that the different 
time and e dependences were confirmed there. 

IV. THEORETICAL ANALYSIS UNDER THE INDEPENDENT CLUSTER 

APPROXIMATION 

Now we turn next to the more powerful analysis starting from a version of the master 
equation, which can provide a full description of the process. This is more easily set up by 
reformulating the process using a column picture, in which a column of height m represents 
a cluster of size m. The map of the problem of particles in a line onto the cluster problem in 
a line of reduced length is shown in Figs. 6a and 6b. The original diffusion processes of free 
and aggregated particles correspond to those shown in Fig. 6b. Since one cluster has two 
edges but corresponds to a single column, the one-particle detachment rate in the column 
picture is 

7 = 2c. (7) 

Finally, in Fig. 7 we illustrate the competition between deposition of a new particle and 
cluster breaking after the formation of a double vacancy in the particle picture, which 
corresponds to the formation of a single vacancy in the cluster picture. 

The deposition process (Fig. 7) leads to the decrease of the total number of clusters and 
of the number of holes between the clusters as time increases. On the other hand, the length 
L of the line in which particles are deposited and diffuse is kept constant. Consequently, 
in order to adopt the column picture, it is necessary to consider that the length L of the 
corresponding column problem decreases in time. These lengths are related as L = L — M, 
where M is the total mass or total number of particles, for periodic boundary conditions. 
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We denote by P t (m) the probability that a randomly chosen cluster (equivalently, col- 
umn) has size m at time t. It is given in terms of the clusters numbers N(m, t) by 

Pt{m) ~ "L^p (8) 
while the length of the lattice in which the column problem is defined varies due to deposition 



as 



L ° (t+ L o( t) L ° (t) = - p ><- 0)l2 - p ><- 0)] - (9) 

Here, -P*(0) is the probability of an empty site in the column problem, which correspond to 
double vacancies in the original particle problem (see Figs. 6 and 7). 

Then the gain/loss from in/out processes provides a master equation which must be 
written in terms of clusters numbers N(m, t) in the general form 

N(m, t + 1)- N{m, t) = L (B m + C m + D m + B m ) , (10) 

where B m comes from the cluster breaking processes, C m comes from the diffusion of aggre- 
gated particles (detachment of particles from clusters), D m comes from the diffusion of free 
particles and B m comes from the deposition processes. The terms in the right-hand side of 
Eq. (10) are thus written in terms of probabilities of cluster masses at time t. 

In an independent interval approximation in which joint probabilities are factorized, they 
are given by 

C m = 7 #(m + 1 - 2)P t (m + 1) + 7 fl(m - l)P t (m - 1) [1 - P t (0) - P t (l)} 

- 7 0(m - 2)P t (m) - 7 P t (m) [1 - P t (0) - P t (l)] , (11) 

An = dB{m - l)P t (m - l)P t (l) + dS m , P t (l) - d6 mtl P t (l) - dP t {m)P t {\) 

= dP t {l) [9{m - l)P t {m - 1) + 5 mfl - 5 m ,i - Pt(m)] , (12) 

B m = P t (0) [29 (m - 2) P t (m - 1) - 29 (m - 1) P t (m) + 5 m>l P t (0) - 25 mfl ] (13) 
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and 

oo 

B m = 29(m - l)P t (0) f( m + r ) p t( m + r ) 



r=l 

-6{m - 2)f(m)P t (m)(m - l)P t (0) - S m>0 P t {0) f> - l)/(p)P t (p) 

oo 

= P t (O)[20(m - 1) 2 /( m + p t(™ + r ) 

r=l 

oo 

-6(m - 2)(m - l)/(m)P t (m) - 5 m , £ (p - l)/(p)P(p)]- (14) 

In Eq. (14), the function f(m) gives the rate for a cluster of mass m to break at each 
internal point, i. e. at each internal connection between two aggregated particles: 

f(m) — k^——^— — k(m - l) a_1 . (15) 
m — 1 

The generating function 

oo 

G t (s) = Y: PtMs m (16) 

m=0 

satisfies 

G t+ i(s)[l - P(0)] 2 - G t (s) =C(B + C + D + B), (17) 
where the operator £ is defined as 

oo 

Cg = £ 9(m)s m . (18) 

m=0 

We thus obtain 

CC = 1 (1 - s) [G t (s) - P t (0) - P t (l)s] + (s - 1) 7 6G t (s), (19) 
s 

with 

6=1-P t (0)-P t (l), (20) 

£L> = dP t (l) (s - 1) [G t (s) - 1] , (21) 
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CB = P t (0) [2 (s - 1) (G t (s) - P t (0)) + sP t (0) - 2] 



(22) 



and 



CB = 2P t (0) J2 



s ™ £ /( m + r)P t (m + r) 



m=l L r=0 

oo 



-P t (0) £ s m (m-l)/(m)P t (m) 



m=2 



-P t (0)£(p-l)/(p)P t (p). 



(23) 



rrt=2 



Eqs. (19), (21) and (22) are the same obtained in Ref. [6] for the model with only deposition 
and diffusion (with d — 1). 

Now defining the operator T so that 



\ OS ) m m 



m), 



(24) 



we may write the contribution of cluster breaking processes to the generating function as 



CB = P(0) 



:i+s) 



s=l 



-P,(0) 



8 = 1 



(25) 



Because deposition slowly fills the system, we expect the configurations to coarsen and 
presumably to go into some scaling asymptotics where mass scales with some power of t, 
and Pt(m) and G t (s) each become one- variable scaling functions. So we look for a long time 
scaling solution of the above equation, in which the finite difference G t+1 (s) — G t (s) in Eq. 
(17) can be taken as a time derivative. The scaling variable will be some combination of t 
(large) and u = 1 — s (small), the latter because large cluster sizes arise from structure in 
Gt(s) at s ~ 1. The variable u is actually conjugate to m. Coarsening will correspond to 
the scale of m as t z , with some power z, in which case the one- variable form will be 



G t {s) =u a F (ut z ), 



(26) 



with some function F. Normalization requires a = and P(0) = 1. In the scaling limit, the 
relationship of the generating function to the probability Pj(m) requires the latter to be of 
the form 
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P tM = \ z 9 (£) , 

with 

F(x) = / g{y)e- x Hy. 
Jo 

Moreover, for s ~ 1 one has sj^ ~ ^ in Eq. (25). 

Thus we obtain the following equation for the one-variable scaling function: 



(27) 



(28) 



(l-^(O)) 1 + 



dt 



- 1 



F(ut z ) 



u 



[F(ut z ) ( a - 1 (l-u + u 2 + ...)) + ( 1 -l)±g(±)+ 1 (l-u + u 2 + ...)±g(0) 

+±g(0) \2uF(ut z ) -(!-«) ±g(0) + MO) - 2 



v 



u=0 



(29) 



where a = 7 - 7 P t (0) + (1 - 7) P t (l). 
The dominant terms in Eq. 29 give 

™F'(x) - 2P t (0)F(x) = l(l) 2 F{x) - 2P t (0) 

+Pt(0)t 



d 



x=0 



x=0. 



+Pt(0)t z 



x 



(30) 



In the scaling limit, the first term at the left hand side (LHS) is of order 1/t. Assuming 
that the density of clusters of zero mass (holes in the column picture) decays as 



(31) 



the second term at the LHS and the second term at the right hand side (RHS) of Eq. (30) 
are both of order 1/t 13 . The first term in the RHS, related to the diffusion of aggregated 
particles, is of order l/t 2z . Since cluster breaking competes with deposition and recalling 
that the coarsening exponent in the problem with only deposition and diffusion is z — 1/2, 
we expect that, in the present system, z < 1/2 for any a. The last term in the RHS, 
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associated with cluster breaking, is of order \/t^~ za \ which clearly dominates over the 
terms which scale as 1/t 13 for any a > 0. Thus, we conclude that the terms associated with 
aggregated particles diffusion (first of the RHS of Eq. 30) and with cluster breaking (last of 
the RHS) are dominant for a > and must cancel each other, while the remaining terms 
are subdominant contribution which must also cancel. From the dominant terms, we obtain 

2z = (3 - za (32) 

and from the subdominant terms we obtain 

= 1, (33) 

which gives the coarsening exponent 

z = 1/ (2 + a) . (34) 

The value of z agrees with that obtained from scaling arguments and in simulations. 
In order to compare the density of single vacancies in the column problem, Pt(0), and the 
density of double vacancies in the original particle problem, p 00 (t) , we consider that the total 
length of the lattice in the column problem, L , is smaller than the length L by a factor of 
the order of the average cluster size (see Figs. 6 and 7), which scales as t z . Consequently 

PooW~^(0)/t 2 ~r (1+2) , (35) 

which also agrees with the scaling arguments of Sec. II and simulation results. The e- or 
7-dependence of average cluster size and densities of vacancies in both pictures, predicted 
in Sees. II and III, also follow from Eq. (30). 

Also notice that, for a < 0, the cluster breaking term of Eq. (30) is not dominant 
anymore and we obtain z = 1/2, as in the problem without that mechanism (Ref. [6]). 

V. DISTRIBUTIONS OF CLUSTER LENGTHS 

An equation for the one- variable scaling function F(x) follows from the balance of the 
dominant terms in Eq. (30), related to diffusion of aggregated particles and cluster breaking: 
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^x 2 F(x) 



+ 



x=0 



x=0. 



(36) 



Physically, it means that the dynamics in long time regime is dominated by diffusion of 
aggregated particles, which forms larger clusters at the expenses of neighboring ones, and 
cluster breaking, which restores a relatively random size distribution after the onset of a 
vacancy in the column picture (double vacancy in the particle picture). Deposition becomes 
a rare process, as expected with the large probabilities of cluster breaking for a > 0, and 
does not contribute to determine the main features of the scaling function. 

Instead of solving Eq. (36) for the function F(x), which would involve fractional deriva- 
tives for non-integer exponents a, we consider the same balance of contributions of diffusion 
of aggregated particles and cluster breaking in the original master equation (10), (11) and 
(14). They can be rewritten in terms of the function g(y) defined in Eq. (27). At this point, 
we must consider that Pt{0) scales as 

(h/k) 



Pt(0) 



t 



(37) 



where ko ~ o(l), since the survival time of those vacancies is inversely proportional to the 
amplitude k of the rate of cluster breaking (Eq. 3). We are thus led to 



7 d 2 g(y) + (k /k) 



t 2z dy 2 



t 



t z kt z(a-l) 



2 I™ dy'y' (a - l) g{y')-y«g{y)-5{y) f° dy'y'^ g(y') 

Jy JO 

= 0. (38) 



For y > 0, Eq. (38) leads to 



W = 2/V + (2 + a) y^g, 



(39) 



with 7! = 7/&0. From this equation we obtain the distribution in the form 



g(y) ~ y Hexp(-Ay v ), 



(40) 



with 



15 



v=l + a/2, fi = 1 - a/4, A ~ 7~ 1/2 , 



(41) 



which suggests that the shapes of the distributions of cluster lengths also depend on the 
cluster breaking exponents. 

This prediction can be compared with simulation results. In all cases, the distributions 
were obtained by taking the length of each cluster only once while spanning the lattice at 
a fixed simulation time. Thus our estimates can be directly compared with the predictions 
from the column picture. On the other hand, if a site average had been done during the 
simulations, a rescaling of the predictions of the column picture would be necessary. 

In Figs. 8a and 8b we show the scaled cluster length distributions for a = 2 and a = 4, 
respectively. From Eqs. (40) and (41), it is expected that In [{et) z P t {m) / decreases 
linearly with y u , with 



([/, — 1/2 and z/ = 2fora = 2,/z = and v = 3 for a = 4). Here a factor e~ z multiplies the 
variable in Eq. (27). 

However, the curved shapes of the plots in Figs. 8a and 8b disagree with the predictions 
of this independent cluster approximation. On the other hand, it was observed that the 
data for all values of a can be fitted by a universal distribution of the form 



This is illustrated in Figs. 9a and 9b for a = 0.5 and a — 1, respectively, in which the 
deviations from Eqs. (40) and (41) are relatively small (indeed, for a = 1, the exponent 
3/2 of the universal distribution is also predicted by Eq. 41), and in Figs. 10a and 10b for 
a = 2 and a — 4, respectively. Least squares fits of the data for the largest lengths (dashed 
lines shifted two units to the right in Figs. 9 and 10) confirm the validity of this universal 
scaling. 

In order to explain this unexpected result, we must pay attention to the details of the 
dynamical process during the coarsening process, which comes from a balance between 



V = m/(et) z 



(42) 




(43) 
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aggregated particle diffusion and cluster breaking, and its connection to diffusion-reaction 
related problems. Turning back to the original particle picture of the problem, during almost 
all the time, the dynamics is equivalent to diffusion of vacancies and scattering of one vacancy 
upon collision of two or more vacancies. This scattering of a vacancy corresponds to the 
breakup of a cluster at a random internal point (see Figs. 2 and 7). Thus, while diffusion of a 
vacancy is the mechanism which leads to the formation of large clusters locally, the collision 
of vacancies restores random distribution of cluster lengths around the average value. 

Focusing on a specific diffusing vacancy, it may be viewed as a particle in an infinite 
sea of traps [12-14], the latter represented by the other vacancies. Subsequent collisions 
approximately correspond to the random walk starting in the origin and the absorption by a 
trap after a certain time - however, contrary to the trapping process, the vacancy is thrown 
away and not absorbed. Thus the probability of a cluster of length x in our model can be 
approximated by the probability of finding a particle at distance x from the origin in the 
trapping problem with an infinite sea of mobile traps (TP). 

In Refs. [13] and [14], the TP was studied analytically and it was shown that the proba- 
bility of the particle surviving at time t, being always located in the interval [—1/2, 1/2] and 
that no trap has entered this region up to time t is given by 



In Eq. (44), p T is the density of traps (of order of the inverse average cluster length, (x) in 
our problem), Ai and A2 are constants of order of a unit and it is assumed that the diffusion 
coefficients of the particles and the traps are both equal to e. This result was used in Ref. 
[13] to estimate a lower bound for the probability of survival of a particle at time t and it 
was shown to converge to the upper bound they also estimated, thus providing an exact 
limit for that probability at long times. 

For a given time t, that distribution is peaked around / = Im, where 




(44) 



A 2 et 



nl/3 



(45) 
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Instead of the broad displacement distribution of a free particle problem, the surviving 
particles at time t are much more probably located around the distance Im from the origin. 
Thus, the probability of finding a surviving particle, which decays in time, may be expressed 
in terms of this most probable distance Im from the origin. In our problem of a diffusing 
vacancy, it corresponds to finding a cluster of length x ~ Im at that time after the last 
collision with another vacancy. This gives a distribution for cluster lengths 

Ai 



P[x) ~ exp 



p T 3/2 x 3/2 - 3p T x 



exp 



C(x/(x)f 2 ]. (46) 



V/ 2 ' 

where we considered that px ~ 1/ (x) and C is a constant. 

The above result agrees with those from simulation for all values of the parameter a. In 
its derivation, a time average was performed, which corresponds to a time average in the 
cluster breaking problem within a large time range for single vacancy diffusion, but small 
for deposition of new particles and coarsening. Indeed, this is the suitable interpretation 
to the fixed-time results of Figs. 8-10, which were obtained from several snapshots of the 
system near the times t given above the plots. 

The discrepancy from the prediction of the theoretical analysis based on the independent 
interval approximation is related to the particular form of the distribution of surviving 
particles displacements in the TP. The former approximation does not capture the correlated 
nature of the rare processes. 

Finally, we observe that the data for several values of a in Figs. 9 and 10, although lying 
approximately in the same region, do not collapse into a single curve, which indicates the 
presence of corrections to the scaling in Eq. (43) which do depend on a. 



VI. SUMMARY AND CONCLUSION 

We studied one-dimensional exclusion models with particle diffusion, reversible attach- 
ment to clusters and deposition mechanisms at large vacancies competing with breakup of 
neighboring clusters. These models aim at representing coarsening of aggregates subject to 
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internal stress, so that the increase of density when there is free space available is restricted 
by breakup that leads to internal relief. 

Different dependencies of the cluster breaking rate on the cluster length were considered, 
involving the cluster breaking exponent a. Simulation results show cluster size growth as x ~ 
t l l( 2+a \ which was explained using heuristic scaling arguments. The analytical treatment 
of the master equation with an independent cluster approximation for joint probabilities 
distributions supports this prediction, as well as the scaling of the density of double vacancies 
obtained numerically. 

Despite the dependence of the coarsening exponent on the exponent a, universal prob- 
ability distributions of large clusters were obtained, in the form P(x) ~ exp (-x 3 / 2 ). This 
result shows limitations of the independent interval approximation for the treatment of rare 
events (cluster lengths much larger than the average value) which emerge from the diffusion 
of vacancies and their survival at long times without being scattered by the collision with 
other vacancies, which corresponds to cluster breaking and restores a random distribution 
of cluster lengths. However, the connection of the problem with the problem of a particle 
diffusing in an infinite sea of mobile particles provides an explanation for that distribu- 
tion, which may be viewed as a snapshot of the system at time intervals during which the 
coarsening process was negligible. 

We expect that the models presented above and the combination of different methods to 
explain their scaling behaviors can be used to understand further non-equilibrium systems. 
A particularly interesting application would be the study of two-dimensional systems subject 
to the same conditions, in order to describe coarsening of adatom islands in cases where there 
is a mismatch of lattice parameters with the substrate and consequent stress or strain of 
those islands. 
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FIG. 1. (a) Diffusion processes of free and aggregated particles at the borders of clusters, with 
the respective diffusion rates, (b) Allowed deposition processes at vacancies with a neighboring 
vacant site, with the respective rate, and the forbidden deposition process, in which the vacancy 
has two occupied neighbors, (c) Three of the possible five configurations after the breakup of the 
cluster at the left side of the double vacancy, with the respective total rate. The breaking positions 
are indicated by a dashed lines. 
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FIG. 2. (a) Three typical neighboring clusters, with lengths of the order of the average cluster 
size x. (b) Growth of cluster A at the expenses of cluster B, which is equivalent to the diffusion of 
the vacancy between them. For simplicity, the effective diffusion of the vacancy between clusters 
B and C was not illustrated, (c) Possible deposition after collision of the two vacancies, extinction 
of cluster B and coarsening of cluster A. 
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FIG. 3. Scaled cluster size as a function of scaled time for the model with a = 1 and aggregated 
particle diffusion rates e = 0.001 (squares) and e = 0.01 (triangles). The solid curve is a parabolic 
fit of the data for e = 0.001. 
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FIG. 4. Scaled cluster size as a function of scaled time for the model with: (a) a = 0.5 and 
e = 0.001; (b) a = 2 and e = 0.001. The solid curves are parabolic fits of the data. 
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FIG. 5. Scaled probability of double vacancies as a function of scaled time for the model with 
a = 1 and aggregated particle diffusion rates e = 0.001 (squares) and e = 0.01 (triangles). The 
solid curve is a parabolic fit of the data for e = 0.001. 
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FIG. 6. (a) Example of particle-hole configuration on a line and the map (dashed arrows) into 
a column problem, (b) The processes of particle detachment from clusters, with rate 7, and of free 
particle diffusion, with rate 1, in the corresponding column picture. Filled squares are the particles 
whose diffusion rates were indicated. 
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(b) 




FIG. 7. (a) Possible events after the formation of a double vacancy in the original particle 
picture: deposition of a new particle (grey circle) and one of the possible events of cluster breakup 
(left cluster) are illustrated, (b) The corresponding events in the column picture. Notice that 
deposition of a new particle (grey square) leads to the decrease of the lattice length in this picture. 
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FIG. 8. Scaled probability of a cluster of size m as a function of the scaled mass according 
to the prediction of the independent interval approximation (Eqs. 40 and 41), for: (a) a = 2, 
e = 0.001, z = 1/4, t = 4 x 10 7 ; (b) a = 4, e = 0.001, z = 1/6, t = 5x 10 8 . 
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FIG. 9. Scaled probability of a cluster of size m as a function of the scaled mass y 3 / 2 , for: (a) 
a = 0.5, e = 0.001, z = 0.4, t = 4 x 10 6 ; (b) a = 1, e = 0.001, z = 1/3, t = 4 x 10 6 . 
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FIG. 10. Scaled probability of a cluster of size m as a function of the scaled mass y 3 / 2 , for: (a) 
a = 2, e = 0.001, z = 1/4, i = 4 x 10 7 ; (b) a = 4, e = 0.001, z = 1/6, t = 5 x 10 8 . 
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